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Abstract — A method for encoding and decoding spectrum 
sliaped binary run-lengtli constrained sequences is described. The 
binary sequences with predefined range of exponential sums are 
introduced. On the base of Cover's enumerative scheme, recur- 
rence relations for calculating the number of these sequences are 
derived. Implementation of encoding and decoding procedures is 
also shown. 

I. Introduction 

Binary sequences with constrained run length of zeros 
are known in Hterature as dk sequences |[T], lIU. In these 
sequences, single ones are separated by at least d, but not 
more then k zeros. A dkr sequence is a dk sequence, ending 
in a run of not more then r trailing zeros. A dklr sequence 
is a dkr sequence, beginning with a run of not more then I 
leading zeros. Among these, sequences with spectral null at 
zero frequency (dc) are notable. This implies a line encoding 
technique using NRZI rules. 

By NRZI encoding [2| we understand mapping the source 
binary sequence x G {0, 1}" to bipolar sequence z £ 
{-1, 1}" such that 

{Zj—\ , Xj — 0, 
— Zj-i, Xj = 1, 

zo = 1. 

The discrete Fourier transform (DFT) of this sequence is 
defined as 

n-l 

i=0 

The usual spectrum shaping requirements are first, a zero- 
mean value for the transmitted data sequence, and second, 
small power content at low frequencies. This second require- 
ment can be written as 

\z*J<M, 

where M is some constant. These zero-mean value sequences 
are called DC-free RLL or DCRLL |[3]. 

Although a study of spectral properties of channel sequences 
was performed by Nyquist H as early as in the late 20s; 
regular papers concerning this problem were appeared in the 
late 60s - early 70s with the contributions coming from 
Gorog Is), Franklin and Pierce [SI, and some other authors. In 
1984 Pierobon Q proved that the finite running digital sum 
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condition is a necessary and sufficient condition for zero mean 
and spectral density vanishing at zero frequency. Later, Marcus 
and Siegel [SJ have expanded the concept of the running digital 
sum to each component of the DFT. This allows extension of 
the spectral null control from zero up to the Nyquist frequency. 

Enumerative encoding and decoding methods for DC-free 
sequences were suggested by Norris and Bloomberg ||9l, 
Immink fTOl, Vasilev lUJ, Braun and Immink [S). 

A method for enumerative encoding these sequences with 
predefined dc component of the DFT Zg = X]j=i ^j, which is 
usually called a digital sum or a charge of the sequence z, was 
suggested in fT2l . Observe that Zq G [— n, n], where Zg admits 
even values whenever n is even and odd values whenever n 
is odd. 

II. Spectrum Shaping 

Let {0, 1}" be the set of all binary sequences of length n 
and let x — (a;i, X2, ■ ■ ■ , Xn) denote a generic element of this 
set. Let .S^in) = {a; e {0, 1}" | satisfies the d, fc, r constraints 
and xi = 1} and let y{n) = {x £ {0, 1}" | satisfies the 
d, k, I, r constraints}. 

Let C,j be the number of sequences from ^(ti); these 
sequences have charge cr = Zq. Using Cover's method ifTsll . 
the number of these sequences can be computed, as is shown 
in lfT2l . using recurrence relation 

r- = /°' n<d + l, 

n \^-W;.'=+l)c-_7^ d+l<n 
with initial conditions 

" 1 0, otherwise. 

We might expand this method to z* . . .z*_]^ components 
of the DFT, but it most likely will not have meaning. In- 
deed, recall that the DFT is a one-to-one transformation. 
If we go from the single spectral component to the vec- 
tor z* = {zq, zl, . . . , z^_i), then for an exact value of 
cr = (dg, (Ji, . . . , (T„_i) S C" we see that C^, can be 
either equal to 1 whenever cr = 2* or equal to otherwise. 
Therefore, we must consider the numbers C^™ for which 
spectral components z*j lie in some area; let this area be a 
ring •s/am,pi m,p2 m ccntcrcd at am with inner radius /9i,,„ e M 
and outer radius p2,,n G K; i.e. the ring £/crm,pi m,p2 m ^ 

■<24„,,pi,„,p2,,„ = {fm : Pl,m < \z*n - CT m\ < P2,m}, ^ 

C. The choice of a ring may be useful in a special case 
when we need only the magnitude range and do not need 



TABLE I 

All Lexicographically Ordered dklr Sequences of 
Length n = 8. Constraints: d = 2, k = 4, 1 = 1, r = 3. 
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By A'^ we denote the lexicographic index of sequence. 




Fig. 1. Paths and spectral components z\ in the complex plane. 

an argument. In this case centre of the ring should coincide 
with the origin of the complex plane and we put cr„, = 0. 

Initially we consider the DFT of n-length sequences. Enu- 
merative scheme implies recurrent calculation of the numbers 
of these sequences. From ^ it follows that each level of 
the recursion diminishes n. Therefore, we must consider two 
different types of lengths; first is the length of the sequence 
and is denoted by 71; second are the lengths of the nested 
subsequences and are denoted by h. Moreover, let the spectral 
properties are given for a large sequence of length n. With 
I and r constraints and the shift theorem, we can obtain this 
ri-length sequence by concatenating ii' G N subsequences of 
length n such that n = Kh. 

Also we intend neither to reject nor even to relax the 
run-length constraints for at least two reasons: first, these 
constraints bound the order of the recurrence relation ^ 
and second, exact values of d and fc constraints should make 
synchronization control easier 

For example. Table |T] shows sequences x, z, spectral com- 
ponents Zq and z\ respectively. The spectral components z\ 
are depicted in Fig. [T] An example of the ring i,pi 1.P2 1 
is superimposed on this image. We also show paths that lead 
from the original to z\. Assuming that spectrum constraints for 
Zq , • ■ • z^ are relaxed, we summarize in Table those se- 
quences for which components z^ lie in the ring i/^s i.pi i,P2 1 ■ 



TABLE II 

All Lexicographically Ordered Spectrum Shaped 
Sequences from Table[i]With ct8 1 = (-2.93, 1.87), 
Pi.i = 1.5,p2,i =2!25. 
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III. The Number of Sequences 

By a;„ = denote a primitive 71th root of unity. Now 
we shall give the following definition. Let < 771 < tt. — 1, 
h < n, z G { — 1,1}, and let 9,^ G C be some constant. Then 
an order j recurrence relation 



t=0 



is called a trigonometric recurrence relation. The constant O^n 
is said to be a linear phase if Om = zuj^"^^. 

This trigonometric recurrence relation is not basically very 
different from the Danielson-Lanczos identity lfT4l . There, an 
77,-point DFT was split into the sum of two DFT; one is 
formed from the even-numbered points, the other from the 
odd-numbered points. Here, we add an exponential sum of 
the prefix to multiplying by the linear phase trigonometric 
recurrence relation of the subsequence. 

Consider bipolar run-length constrained sequences of length 
n and having an 771th trigonometric recurrence relation ring 
-^Cfi m,pi m,p2 n, ccntcred at crf^„i with inner radius pi^m 
and outer radius p2,m', here the ring -s^crs m,pi m,P2 m a 

set ^o-s,„,Pi,„..p2,„x = {0"n,m : Pl,m < ~ CTfi.ml < 

P2,m}, <„ e C. 

We compute the number of n-length subsequences of z. 
Let (7^"j"'''i."'''2,n, jj^g number of these subsequences, 
which begin with one. Let (^^"^"'''i.'"'''^,^ number of 

these sequences, which begin with a leading run of zeros. 

Since an internal run of zeros succeeds a leading 
run of zeros, we see that the leading constraint I does 

For convenience, below, under 
{d,k, r) and under 



not affect C^X^ 



,P2,n 



'Pi, 



,Pl,m,P2/ 



,Pl,»7.,P2.' 



Proposition 

^0'S_m,pi^„,,p; 

'ii ft fn 



we imply C^y;^ 
we similarly imply C^"^"l 

1. The numbers C° 
can be obtained as: 



,Pl,m,P2,„ 



and 



-\min(n,/i; + l) ^i<7n-j.'tn:Pl,'tn:P2,'. 

jj—d-\-l n.h—j^ni 



n < d+l, 
bfi,m, d+\ <n. 



where 



-m(j-l)/2: 



(2) 



(3) 



initial condition 

_ j 1, n < r + 1 and pi^,n < \crn,7n - o-fi,,„| < P2,, 
I U, otherwise, 

and additional condition 

-1, r + l<h<k + l 

and pi^rn < ICfi^m - < P2,i 



where 



O'fi.m — —LO. 



If a leading series is running, then 

min(n,/) 

n,h,Tn / ^ n,n— J 



(4) 



1, k + 1 < n <r + 1 (5) 

and pi^rn < \<7h,m - ^n,m\ < P2,m, 

0, Otherwise, 



(6) 



Multiplying (fTTT i by —1, we obtain (|6]i. 

Now we must prove the recurrence relation in Accord- 
ing to Cover's enumerative method ifTSll . we build the recursion 
by the following way. Let us consider a possible run of zeros, 
which follows the leading one, as a prefix for the following 
subsequences beginning also with one. Assuming the length j 
of the prefix grows from d + 1 to min(7i, fc + 1) and weight 
of this prefix equals one, we can consider a concatenation of 
the prefix and the following subsequences as 



d+l<j<k+l 



1 0. 



1 0. 



-1 -1 . . . -1 



In fact, if the subsequences of length n begin with one, then 
the mth trigonometric recurrence relation for cr ,~i,„i can be 
obtained as 



-m0-l)/2: 



Sin 



additional condition 

- J"' 

I 0, Otherwise 



r < h < I and 



(7) 



(8) 



Here d > 0, k > d, I > 0, r > 0. 

Proof: First, consider initial and additional conditions (|4]i 
and (|5]l. In the case of n < r + 1, there is only a trailing run 
in the subsequence. This trailing subsequence of x consists 
of one and n — 1 zeros. The NRZI rule takes each term of 
the trailing subsequence to —1. It defines the mth exponential 
sum (Tfi_,„, which corresponds to this subsequence, as follows: 




Therefore, it gives us the only allowed sequence which length 
lies in the interval [l,r + 1] and an^m lies in the ring 

'^(^fi m-Pl m-P2 m' 



Pl,m < \crfi,r. 



(9) 



Consider a finite geometric series; then recall the derivation 
of the Dirichlet kernel trigonometric identity and obtain 



n-l 
t=0 



1 - w„ 



1 - Wn' 



-m(n-l)/2^ 



sm — TT 



(10) 
(11) 



where the first term is the riith exponential sum of the 
prefix and (Tfj_j is the rnth trigonometric recurrence relation 
corresponding to following subsequence; this subsequence also 
begins with one. The phase factor 6m can be found by the shift 
theorem as follows: Before concatenating a prefix, for the first 
term zi of z, we have 

Since there is just a one in the prefix, it follows that this prefix 
changes a sign of each term in the following subsequence. 
Therefore, after concatenating a prefix of length j, for the 
same term we obtain 

Since zi corresponds to Xi and Zj+i corresponds to Xj^i, 
where xi and Xj-\.i are the same term, it follows that 



Therefore 



J-1 

E —mt —mj 



(12) 



Consider a finite geometric series in this equation; then 

1 - c^-™^' 



Vc.-"* = 



-m(j-l)/2: 



and substituting it for the sum in (fT2] l. we get 



-a~ _a;-™(J-i)/2 



sm 



Multiplying both sides by w™-', we obtain 

The other case is when a leading series of zeros is running. 
First, we also consider the additional condition. In the case 
of n < min(^, r), the leading run of zeros is the trailing one. 
This also gives us the only allowed sequence which length lies 
in the interval [0, min(^, r)]. 

0<n<min(/,r) 



We state without proof the following proposition. 



-1 -mt 



Using ( [Tol l we get 



1 - 0J~ 



i — OJn 

^.-„Uf.-i)/2 ^in(if7r) 
sin(^.)- 

Substituting it for crs m in (|9]l, we obtain the condition in (O. 

In the case of nonzero weight, there exist only zero weight 
prefixes which length lies in the interval [0,^]. Subsequences 
beginning with one follow the prefixes; therefore, we can 
consider a concatenation of this prefix and the following 
subsequences as 



0<j<l 




Thus, the mth expression for an,m is 

t=0 

Recall the reasoning that led us to (|3]l. Arguing as there, and 
taking into account that the prefix now consists of all zeros, 

we get dm ~ "^n ™d 



-co-"'' a. 



I- iN/oSin(— tt) 

- m ( J — 1 ) / 2 \ n / 



Multiplying both sides by — w™^, we obtain (|7|i- □ 
Now consider a set Aa-^^.p^^p^ of the rings 

Let (Tfj = (crfi,o,o-na, • ■ • e C" be a vector of 

centres of these rings, pi ~ {pi.o, Pi,i, • • ■ , Pi.n^i) G and 
Pi = {P2,Q, P2,i, ■ ■ ■ , P2,n-i) G R" vectors of theirs inner 
and outer radii. 

Then consider those n-length subsequences from S'{n) 
whose vectors of the exponential sums belong to Ao-j^.p^.p^. 

Let C^};^^-^^ and C^};^^'^^ be the number of these 
sequences, which begin with one and with a leading run of 
zeros respectively. 

As above, under C^"~^^'^^ we imply C^"^^^'^^ {d, fc, r) and 



Proposition 2. The numbers C,'^'^''^^'''^ and C'^"f{^^'^^ 
be expressed as: 



c: 



n < d+1, 



(13) 

where components (Jn-j,m of the vector (Tfi-j ~ 
{an-j,o, an-j^i, . . . , crfj_j.„_i) are defined by (|3]l. Here initial 
condition 

{1, n < 7' + 1 and 6 is true, 
(14) 
0, otherwise, 

and additional condition 

— 1, r + l<7i<fc+l and 6 is true, 

bf, ~ ( 1, k + l<h<r + l and 6 is true, (15) 

0, otherwise, 

where the indicator function 5 is the logical conjunction of n 
statements each of which predicates of the {'m)th exponential 
sum of trailing run (|6]l lies in the ring ^ ^ .p^ i.e. 

71-1 

^ ~ /\ {pi,m < \0'A,m - <^h,m\ < P2,', 
rn—0 



If a leading series is running, then 

c: 



min{n,l) 

where components ^fi-j,m of the vector fTn^j ~ 
{^n-j.Q,crn-j,i, ■ ■ ■ ,(^n-j.n-i) <^f£ defined by O. Here ad- 
ditional condition 



bn = 



^1, r < n < I and S is true, 
0, otherwise. 



under C, 



we similarly imply C„ V^^''^^ {d,k,l,r) 



where the indicator function d is the logical conjunction of n 
statements each of which predicates of the {m)th exponential 
sum of trailing run —an,m H^s in the ring s^a^ rn.pi m.pi m' 
i.e. 

n-l 

^ = /\ (^Pl,m < \crA,m + <^h,m\ < P2,', 

IV. Algorithms for Encoding and Decoding 
Spectrum Shaped dklr Sequences 

Now let denotes a set of spectrum shaped binary run- 
length constrained sequences x ~ {xi, X2, ■ . ■ , Xn) of length 
n. Let the set ^ be ordered lexicographically. From it follows 
that the lexicographic index N{x) G {N, 0} of a; G is given 
by 

fl 

N{x)^Y.'^,W{p), 

where W{p) denotes the number of sequences in 5^ with 
given prefix p = {xi,X2, ■ . ■ ,Xj_i,0). 



The decoding algorithm, for given sequence x, find its 
lexicographic index N{x) < This is done by successive 
approximation method using W{p) as the weight of term Xj. 

By aj (p) denote the number of trailing zeros of the prefix 
p. By Vj-i = X^ti denote the weight of this prefix. 
Since p is the prefix of x, it follows that subsequence x = 
(ccj, ccj+i, . . . , Xn) is the rest of x, and Ij is the leading run 
of zeros in this subsequence. We define Ij as the complement 
of aj (p) in I (for a leading run of zeros) or in k (for the other 
runs of zeros) as follows; 



k~aj{p), otherwise. 



Now we can compute the number of the spectrum shaped 
sequences Wa-f, (p) as 



W{p) 



71, n— J 



0, 



otherwise. 



By hfi we take into account a trailing run as follows: 

1, r — aj{p) < n — j < min(Zj, r) and S is true, 



bn 



0, otherwise. 



where the indicator function S is the logical conjunction of n 
statements each of which predicates of the mth exponential 
sum of the trailing run 



i(n-i-l)/ 



— —n 

n J 



lies in the ring ^(T^._j,„,pi,„,p2,„. i-S- 



^ = /\ (pi,m < Icrn-j.m — d'a-j.ml < 



m=0 



where components (Jn-j.in of the vector cra-j ~ 
(cTh-j^o, <^A-],i, • • ■ , crfi-j>i-i) are defined as 

We introduce the next variables: a, w, c = (co,ci,..., 
Cn-i), which correspond to aj{p), Vj-i, crj-i- 

N:=0; a := 1; w := 0; c 0; 
for j := 1 to 71 do 
Get W{p); 
if = 1 then 

N:^N + W{p)] a := 1; w := u; + 1; 
else 

a ;= a + 1; 
end if . . . else 
for 771 to n — 1 do 

c™ c„„ + (-1)"W"^^-'^ 
end for 
end for. 

The encoding (inverse) algorithm, for given lexicographic 
index N{x) < \5^\, find the corresponding x. 



a := 1; w := 0; c := 0; 
for j := 1 to n do 

Get W{p); 

it N > Wip) then 

N:=N-W{p); xj := 1; a := 1; 11}:= w + 1; 

else 

Xj := 0; a := a + 1; 
end if . . . else 
for 771 := to 71 — 1 do 

Cm '•— Cjji + ( 1) LOn j 

end for 
end for. 

V. Conclusion 

We presented an enumerative approach for constructing 
binary run-length limited sequences to satisfy spectral con- 
straints. By considering the spectral components of the DFT 
in the complex plane, we got the possibility to use Cover's 
enumerative scheme for encoding these spectrum shaped se- 
quences. First, we defined the sequences whose 7nth spectral 
components lie in a certain ring. Secondly, we proposed recur- 
rence relations for calculating the number of such sequences. 
Then, we expanded our method for vectors of all spectral 
components. Finally, we suggested algorithms for encoding 
and decoding these sequences. 
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